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Abstract 

A phenomenological theory is presented for two-dimensional quantum liquids 
in terms of the Fermi surface geometry. It is shown that there is a one-to-one 
correspondence between the properties of an interacting electron system and 
its corresponding Fermi surface. By doing this, the concept of Fermi surface is 
generalized to include different topologies. It is shown that for a Fermi liquid 
the corresponding Fermi surface is rough. In the presence of a condensate, 
the Fermi surface is faceted, while for a ferromagnetc instability, the Fermi 
surface becomes a frozen solid. I also determine the surface tension, the step 
free energy, low lying excitations and other surface and transport properties 
of the Fermi surface. The different transitions between these phases are also 
determined. A non-Fermi liquid phase is shown to be a pre-roughening state 
of the Fermi surface, the properties of which are briefly analyzed. 
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I. INTRODUCTION 



I present a completely new semi-classical approach to the quantum liquid theory in terms 
of true surface properties and geometry of the Fermi surface. Within such a description, it 
can be shown that there is a one-to-one correspondence between the properties of a two- 
dimensional interacting electron system and its one-dimensional Fermi surface. Contrary to 
popular belief and expectation, I will also prove that the Fermi liquid theory corresponds 
to a rough and not smooth Fermi surface. Some of the evidence can be applied to a three- 
dimensional interacting electron system. 

The bulk of the paper is dedicated to the detailed analysis of the rough phase and the 
corresponding roughening transition. For completeness, I present in Sections [II] and |T| the 
general theory of equilibrium surfaces and the roughening transition. In addition, the theory 
of the more particular Pokrovsky-Talapov transition is presented in the Appendix. This 
standard surface physics theory is applied to Fermi surfaces in Section |V[ while in Section 
[V| the possible transition between a rough and a flat faceted Fermi surface is described. Due 
to the present interest in non-Fermi liquid theories, I will also give a possible scenario in 
Section [Vl| when and how such a phase could appear. 

II. EQUILIBRIUM SURFACES 

The theory of the equilibrium surfaces dates back to the celebrated work of Wulff (Wulff 
1901) and by the time of Herring's 1951 paper (Herring 1951) the theory was truly under- 
stood. The equilibrium surface shape is the shape that minimizes the total free energy of 
the surface at a fixed volume. This constrained minimization was originally solved by Wulff, 
which is an elegant construction and valued textbook material. To find the equilibrium 
surface shape, follow the rules (Wulff 1901): i) make a polar plot (the so-called Wulff plot) 
of the surface tension a(n) as a function of n, ii) starting at the origin, draw the radius in 
direction n out to the Wulff plot, then Hi) construct planes perpendicular to n's through the 
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intersection of the radius with the Wulff plot. The interior envelope of this family of planes 
is always a convex figure whose shape is that of the equilibrium surface^. 

In the following, I present the analytic derivation of the equilibrium surface shape based 
on the Landau and Andreev derivation (Landau and Lifshitz 1980), which is not the most 
general (it relies on a(n) being piecewise differentiable) but it has the advantage of being 
clear and easy to follow. 

As noted previously, the Wulff construction always produces a convex shape, which can 
be divided into a top and its corresponding mirror-image bottom surface, with respect to the 
+z and —z direction. For simplicity, hereafter I only deal with the top surfaces, as all the 
obtained results can be mapped to the bottom one. Through out this paper I will use the 
notations: z(x,y) is the height of the surface with slopes dz/dx = h(x) and dz/dy = h(y), 
steps on the surface will have height H and width L with corresponding tilting angle 9. 
With these notations the total free energy of the surface is: / dxdyf(n), with 



f(h) = a(h) yjl + h 2 (x) + h 2 (y) (1) 

being the surface free energy per unit area. Minimizing the total surface free energy keeping 
the volume fixed is identical (Landau and Lifshitz 1980) to solving the variational equation 

5 J dxdy [f - 2\z] = , (2) 

where A is a Lagrange multiplier. Eq. @ transforms into 

[ d d \ ^ ^ Z + ^ ^ Z 2X5 \ (3) 
J \-dh{x) dx dh(y) dy -I 

By assuming its variation to be equal to zero, for any height variation 5z, Eq. (Q) can be 

partially integrated. The obtained result, again, must be valid for any 5z, accordingly a 

simple integration with a proper choice of the integration constant gives: 

d df d df 2A 
dx dh(x) dy dh(y) 



1 The overall magnitude has no significance. 
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The solution of Eq. (|j) is (Landau and Lifshitz 1980): 

/ = A [z - xh(x) - yh(y)} . (5) 

Eq. (|5|) is non other than the Wulff plane perpendicular to n and the interface profile will 
have the form: 

I Of 1 df 1 r df df 

x = — — , y = — — and z = — J — n(x) — h{y) 



X dh(x) ' * X dh(y) ' X L v 7 dh(x) Ka/ dh(y) 

(6) 

The shape of the surface z(x, y) is obtained as a Legendre transform of the surface tension 
per unit area, and the surface free energy is the Legendre transform of the surface shape: 

Oz Oz 

f — z — x — — y — . (7) 

dx dy 

Note, that the above construction is so general that it applies to any equilibrium interface. 
As an example I present in Fig. 1 the equilibrium shape of a square with its Wulff plot. It is 
important to emphasize that in first approximation one would consider the surface tension 
of a square to be vanishing, as the surface radius is infinitely large. However, this is not 
the case. For a proper evaluation of the surface tension we must count all the broken bonds 
at the surface, or simply the number of bonds to which a periodic boundary condition is 
applied (number of periodic bonds). Consider, for simplicity, that for the square in Fig. 1 
the interactions in z and x directions are equal to unity, then the calculated surface tension 
is: a in) = I sin 01 + 1 cos 9 1 . 



III. THE ROUGHENING TRANSITION 



The disappearance of a cusp, sharp edges, canonical points, etc from the Wulff plot are all 
related to phase transitions. From these, the most studied is connected to the disappearance 
of cusps. Since the first theoretical attempts by Burton and Cabrera in 1949 (Burton 
and Cabrera 1949), this singularity was usually called the roughening transition or the 
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phenomena of surface melting. These two denominations refer to the same thermodynamic 
transition. Historically roughening transition is usually used in the theoretical works where 
the emphasis is on the critical phenomena aspect of the phases transition, while surface 
melting is an approach from the point of view of surface physics and experimentalists. This 
may generate some confusion, but everything becomes clear if we follow the argument of 
Burton et ed. (Burton and Cabrera 1949; Burton, Cabrera and Frank 1951) outlined in 
the following: the static properties of a two-dimensional surface are well described by a 
Hamiltonian0 

H = \Hi - Hj\ , (8) 

where the sum runs over all pairs of nearest-neighbour sites. With the simplest choice for 
the height variables, Hi = {0,1}, Eq. fl8|) is identical to the two-dimensional Ising model 
(S{ = 2Hi — 1) which undergoes a second-order phase transition at k B T c = 2.27 J , below 
which (Hi) < 1/2 and above which (Hi) = 1/2. For the surface this means that above 
T c the surface is rough, as the surface would consist of 50% each atoms and holes, but the 
nucleation barrier for two-dimensional nucleation would vanish, accordingly the surface is 
melted. 

The presence of a cusp in the Wulff plot, is due to the appearance of an absolute value in 
the expansion of the surface free energy or surface tension on the smooth side of a surface, 
e.g., for a one-dimensional interface: 

f[h(x)] = /(0) + \h(x)\f B + ... , (9) 

or 

a(6) = a + \6\ a% + (10) 



2 Some authors use the disordering transition denomination. 

3 This is the well-known Solid-on-Solid model for a crystal- vapor interface, or the "Kossel" -crystal. 
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respectively. The coefficient / s of called the step free energy is the energy cost of 

creating a unit length of step or ledge of unit height. When is finite, adjacent steps are 

well separated and the energy per unit step length is well defined. The size of a horizontal 
facet is proportional to the jump df/dh(x), or equivalently, to the jump in da/ 36 [as 
tan# = h(x)} at the cusp which follows immediately from Eqs. (^) and (Q): when h(x) crosses 
zero, x jumps from +f s /\ to —f s /X, indicating that the facet size is 2/ s /A. Accordingly, the 
roughening transition is characterized macroscopically by the disappearance of a cusp from 
the Wulff plot or governed microscopically by the vanishing of the step free energy. The 
Wulff plot will be cusped for facets with smooth surface, while the cusp disappears from the 
Wulff plot of rough rounded surfaces. 

Usually, the singularity at this transition is in the same universality class as the 
Kosterlitz-Thouless (1973) transition. The smooth interfacial phase maps to the high- 
temperature, exponential XY phase, while the rough interfacial phase maps to the low- 
temperature algebraic (Gaussian) XY phase. Or in sine-Gordon language, the rough phase 
maps into a non-interacting, or Tomonaga-Luttinger model, while the smooth surface into 
a Luther-Emery model. However, there are other types of continuous transitions, which can 
give rise to a rough surface, such as Ising, Potts, Pokrovsky-Talapov (1978) or chiral that are 
more frequent, depending upon symmetry, number of states, interaction parameters, fields, 
chirality, etc. An analysis of the critical behaviour of these transitions is not the purpose of 
the present work. Such a microscopic approach to the rough surfaces in general and to the 
Fermi surface in particular will be published elsewhere (Gulacsi 1996). 

Most real crystals grow under conditions that are far from equilibrium. The application 
of the present theory for a real surface is limited by the requirement of equilibrium. From 
this point of view the Fermi surface is an ideal surface as by construction^ it is in equilibrium. 

4 The fixed volume condition is guaranteed by the Luttinger theorem. In most of the cases through 
out this work, the Luttinger theorem is local a la Haldane (Haldane 1994). 
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IV. SURFACE PROPERTIES OF FERMI LIQUIDS 



The Fermi liquid picture is both simple and profound. With a few basic assumptions 
it describes accurately the low-energy, long-wave length behaviour of interacting systems in 
particular and the thermodynamic properties of the whole system in general. The theory 
assumes that there is a one-to-one correspondence between the low-energy excitations of a 
non-interacting system and the low-energy excitations of the interacting system while pre- 
serving the quantum numbers^} That is, if the non-interacting system is characterized by 
some distribution function n(k) of the bare particles, by switching on the interaction adia- 
batically the eigenstates of the interacting system, called quasiparticles, are also described 
by same function n(k). Accordingly, the excitation of the system is simply measured by the 
departure 5n(k) from the ground-state distribution: 

5n{k) = n(k) - n°(k) . (11) 

For the ideal system, there exists a simple linear relation between the energy of a given 
state and the corresponding distribution function. When particle interaction is taken into 
account, the relation between the state energy, E, and the quasiparticle distribution function, 
n(k), becomes much more complicated. It may be expressed in a functional form, £7[n(k)], 
which one cannot in general specify explicitly. If, however, n(k) is sufficiently close to the 
ground-state distribution n°(k), we can carry out a Taylor expansion of this functional. On 
writing n(k) in the form Eq. (|TT|) and taking 5n(k) to be small, or to extend over a small 
region in momentum space, the Landau free energy (Landau 1957) becomes: 

F - F = ]T[e(k) - fi] 5n(k) + W /(k, k') 5n(k) 5n(k') + ... , (12) 

k 1 k,k' 

which is the heart of the phenomenological theory of Fermi liquids (Pines and Nozieres 1988) 
with /(k, k') representing the effective interaction between quasiparticles. 



5 For simplicity in the following and through out the present work only the spinless case will be 
discussed. 
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It is not my purpose to further analyze Eq. ([12]) from the Fermi liquid point of view[], 
but only to establish its surface properties as a semi-classical interface. For simplicity, I 
restrict ourselves to one- dimensional Fermi surfaces[] (two-dimensional systems). It should 
be noted that some formulas or demonstrations that appear in the text are also valid or 
easily generalized for a two-dimensional Fermi surface (three-dimensional systems). The 
reader will be prompted when these examples arise. 

Based on the definitions presented in Sections |T| and [ill] it is immediate that Eq. (|i~2"D 



describes a rough surface because the step free energy is zercf]. As noted already by Pines and 
Nozieres (1988) Eq. ([12]) looks like a Taylor expansion, but it is not. With the notations from 
the previous sections and only if \i and 5n(k) are finite, F — F defined in Eq. (|T2| ) is of the 
order /i 2 (fca;)(] because e(k) — /i oc h{k x ) and 5n(k) oc h{k x ). Accordingly, the corresponding 
step free energy, / s , is strictly zero and by definition the surface is rough. The roughening 
transition is defined by the vanishing of the step free energy, and any smooth surface must 
characterized by Eq. (|9|). The only case where such a simple proof would not work is the 
case /i = 0, which would correspond to a half-filled band systems, a case which I will analyze 
later. 

This proof should work equally well for two-dimensional Fermi surfaces for the non- 



The excitation of Eq. (|I^), i.e. the sound modes and the particle-hole continuum, will be briefly 



analyzed in Section VII 



7 The surface under study, being defined in the momentum space, the local curvilinear surface 
coordinate is k x , k z (k x ) denotes the local height of the surface with slope dk z /dk x = h{k x ). The 
tilting angle remains 9. It is also important to emphasize that the steps are associated with the 
world lines of fermions (Pokrovsky-Talapov 1978). 



By construction Eq. (12) describes density fluctuations of a liquid surface, rather than elastic 



waves propagating on a smooth facet. 
9 See, Eq. ( |i~9| ) for an expansion of a rough surface. 
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interacting case. However, some problems will appear if we consider the two-dimensional 
Fermi surface, a quantum interface case, which is analyzed in detail by Gulacsi (1996)0. For 
a one-dimensional Fermi surface, the above argument is based on a nearest-neighbour type 
effective interaction. However, at large length scales the discreteness of the step heights is 
irrelevant, the surface always can be described by a Gaussian model, having the properties of 
a transverse vibrating elastic network. Due to this, the obtained transition (for details see, 
Section |V|) is not destroyed if the interaction between quasiparticles is of short-range type 
or even has the form /(k, k') ~ g/\k x — k'J^ with 7 < 2 and g > 0. From these interactions 
the g/\k x — k' x \ 2 form is the most probable one (Ashcroft and Mermin 1976) based on a 
Hartree-Fock theory, or simply because the elastic or bipolar interactions between steps on 
an equilibrium surface have this form (see, Section [V]). 

The Fermi surface is indeed rough for the non-interacting case, and it can be directly 
proven by the Wulff plot. Equations ([J) through Eq. (0) enable one to reconstruct the Wulff 
plot from the equilibrium surface shape. For an isotropic and non-interacting system with 
ep = ^ 2 |kir| 2 /2m we immediately obtain that 01(6) = ep becomes independent of 9, as in the 
case for an ordinary fluid. Accordingly it has no cusps; that is, the Fermi surface is rough. 
In this simple case, following Nozieres (1991), we can determine the Lagrange multiplier of 
the Fermi surface appearing in Eqs. ©: A = v F /4m. 

On a lattice^ the half-filled and doped case must be treated separately. Namely, for 

10 Fisher and Weeks (1983) and Fradkin (1983) and others have contested the Andreev and Pashkin 
(1978) argument about the roughness of the quantum liquid surfaces at zero temperature. However, 
new arguments appear in the favor of the roughness, namely as presented in Footnote 14, the width 
of the steps on a quantum surface at T — > are very large, so the quantum effects are becoming 
smaller, accordingly the quantum surface becomes classical before the complete derealization of 
the interface. 

n For simplicity consider a square lattice. 
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H = — tH(i,j)( c l c j + h- c -) + l^Hi c t c i, with [i = —2t5, where 5 denotes the hole doping, using 
Eq. ®, the Wulff plot corresponding to the Fermi surface /J = 2tJ2i cos ki will be smooth 
for /i 7^ and cusped for \x = 0. A similar argument used in the previous paragraphs is 
that Sn(k) oc h(k x ) and also e(k) — /x oc h(k x ) and for /i ^ 0, i.e., F — Fq oc h 2 (k x ), while 
for /j = only 5n(k) oc h(k x ), so F — F oc h(k x ), accordingly the step free energy is finite. 
In this case, the Fermi surface and its Wulff plot are identical to Fig. 1 (rotated by tt/2). 
That is, the Fermi surface is faceted and the facets are smooth. However, in this case, 
the appearance of the facets can be attributed to the presence of the van Hove singularity. 
The step free energy is finite but small close to the van Hove singularity, increasing as we 
approach the singularity. This behaviour is caused by the vanishing group velocity, |Ve(k)|, 
which explicitly shows up in evaluating Eq. (|7|). The step free energy, however, cannot grow 
infinitely large, is bounded by the crystal structure, but produces a facet with the largest 
possible flat surface, and this is the Fermi surface of the half-filled band case. The above 
arguments for the half-filled band case are also valid for a two-dimensional Fermi surface. 

In conclusion, in a non-interacting system the change of the topology of the Fermi surface, 
i.e. the appearance of flat facets can be attributed to the presence of the van Hove singularity. 
For an interacting case it is difficult to determine the Wulff plot because we do not know 
the exact form of the self energy, £R e (k, u). However, it can be shown that the Fermi 
liquid picture in general, and Eq. ( O ) in particular, is not valid close to half filling. For 
example, for the two-dimensional Hubbard model, second-order perturbation theory gives 
for the /(k, k') function (Honner 1995): 



lA AU 2 Ak , 4 sin Ak + Vl6 sin 2 Ak - A 2 e . . 

/(k,k) = , = In , , (13) 

Ntn Vl6 sin 2 Ak - A 2 e 4 sin A A; - Vl6 sin 2 A A; - A 2 e 

where 2Ak = k x — k' x = k y — k' y and Ae = cos k' x — cos k x + cos k' y — cos k y . The /(k, k') from 
Eq. (|i~3| ) is bounded only in the vicinity of (tt/2, tt/2) with a maximum value of 2U 2 /Nt7r 2 
and is divergent as we approach the van Hove singularity. That is, the adiabatic switching 
on the interaction picture of the Fermi liquid theory is no longer valid. 

For an interacting electron system, the half-filled band case is mostly due to perfect 
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nesting, in the presence of which the Fermi surface becomes unstable to the formation of 
long-range order (e.g., charge- density- wave or spin-density- wave), or Mott insulator, etc., 
due to which a charge or spin gap (depending on the interactions) will open. This causes 
the divergence of the /(k, k') function. In this case, the Fermi surface is not defined in the 
sense of Eq. (|TT|) or flT2p. Within the present approach, based on the Fermi surface geometry, 
we can generalize its definition with the understanding that only its topology changes, which 
means that the Fermi surface becomes faceted. Previously, the Fermi surface was rounded, 
while now it has flat, faceted regions. The presence of a facet would immediately indicate 
that the Fermi surface cannot be excited by any small amount of energy. A minimum energy 
is required (the step free energy defined in Section [□]) to create steps on a flat and smooth 
facet and this minimum energy signals the presence of a gap in the excitation spectrum. 
In the next Section the width of the facet is determined, in first approximation, to be 
proportional to A/vf, where A is a mean-field gap and Vf is the non-interacting Fermi 
velocity. Accordingly, the appearance of these facets on the Fermi surface is equivalent to 
an instability of the Fermi surface by which a condensate will emerge in the interacting 
electron system. The presence of facets in a two-dimensional Fermi surface will equally 
mean an instability of the Fermi surface towards a quasiparticle gap opening. 



V. THE FACETING TRANSITION 

The half-filled band analyzed in the previous Section is a particular case where all the 
Fermi surface is flat. Generally, facets appear on a surface gradually, meaning their width 
increases with variation of some parameters (doping, external magnetic or electric field) or 
temperature. In any case, the appearance of facets on the Fermi surface will always lead to 
nesting (due to the symmetries of the Brillouin zone) usually a partial one. As the surface 
of the facets increases, nesting becomes stronger and stronger until it reaches a completely 
nested Fermi surface similar to the half-filled band case. The transition from the rough 
surface into a faceted one, is a roughening transition by which the topology of the Fermi 
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surface changes completely, and the properties of which are determined in the following. 

The surface free energy of a facet is given in Eq. ([]). However, this term represents only 
non-interacting steps. To describe properly the faceting transition we need the knowledge 
of the interactions between the steps present on the surface. As the present discussion is 
a phenomenological one, we do not intend to subtract these interactions from perturbation 
theories but rather to analyze them from a surface physics point of view. There are several 
main sources of interaction on a surface. The first is the fact that we do not allow steps to 
cross, or create overhangs, as these states would require higher energy. Meaning, the steps 
will be restricted resulting in an entropic interaction. The contribution of such an entropic 
repulsion, we can understand by following Gruber and Mulling (1967): let us suppose for 
a moment that a single step is constrained to move between two walls situated at k x = ±&£. 
The presence of the walls produces an entropic repulsion between the wall and step. The 
probability density P(n) for the step passing true k x = n is (Gruber and Mullins 1967): 

p(n > = ir os si • (14) 

and the free energy of the step becomes^: / s = / s (0) + const jk\. Now, replacing the fixed 
boundaries of the walls by neighbouring steps results in: / s = / s (0) + const h 2 (k x ), which 
introduced in Eq. (|S|) will give an ^(A^)] h 2 (k x ) contribution. 

Short-range interactions (Jayaprakash, Rottman and Saam 1984) cannot alter this be- 
haviour. Namely, in a mean-field theory the contribution of a short-range interaction V{k x ) 
will be {1/L) J dk x V{k x )h{k x )G{k x ) , where G(k x ) is the non-interacting pair correlation 
function, G{k x ) = h(k x ){l — sm 2 [iik x h{k x )/ L]/[7ik x h{k x )/ L] 2 }. If V is smaller than the dis- 
tance between steps or the step widthF], G(k x ) may be expanded and the contribution of a 



12 The original proof was for a terrace-step-kink (TSK) type step. 

13 If P{n) would be a constant, than the step staying away from the walls would decrease its 
entropy and thereby its free energy with 1/L. 
14 The step width in the T — » limit can be related (Balibar, Guthmann and Rolley 1993) to the 
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short-range interaction will be of order h A (k x ). 

Other interactions (see e.g., Beijeren, van and Nolden 1987), e.g., elastic or dipole in- 
teractions, assume the g/\k x — k' x \ 2 form on a surface with free energy contribution in the 
continuum limit with P(k x ) given in Eq. (0), true the whole length of a facet at small h{k x ) 
is: 

~ / dk x - * ' r : cos 2 nk x h{k x ) ~ \h(k x )\h 2 (k x ) . (15) 

J-(2\h(k x )\)-i k x + l/2\h{k x )\ x K xJ K J 

Accordingly, on a facet the surface free energy from Eq. (|T^ ) will have the form: 



f[h(k x )} = /(0) + f B \h(k x )\ + f iQt .\h(k x )\h 2 (k x ) + ... , (16) 

where in f int we can accumulate the effects of all the interactions present on the surface. 

At this point, it is important to emphasize the following: i) All type of interactions on a 
surface are giving a contribution of h x in the expression of the surface free energy. Naively, 
one might expect the free energy due to interactions to be of the h 2 (k x ) order, however it is 
not. A contribution of the h 2 (k x ) type occurs only in a rough surface, see Eq. (|19|) or Eq. 
(|T2|) and the discussions following it. ii) The presence of the interaction term in Eq. ( |T6"D 
is crucial for the faceting transition. The f int = corresponds to a mean-field type critical 
bahaviour (see e.g., Rottman and Wortis 1984) which may containe a tricritical pointf^. 

In the following I determine the interface profile. Continuing the arguments from Section 
PI , when h(k x ) crosses zero, the horizontal width of the facet, k x , jumps from +f a /X to 
—f s /X, that is the facet width will be proportional to f s /X. In a mean-field approach the 



spectrum has the ye 2 (k) + A 2 (k) form. Using the continuum result from Section [TV] for the 
Lagrange multiplier: A = vp/4, in first approximation it is obtained that the Fermi surface 



so-called pinning or coupling strength of an effective sine-Gordon model, accordingly the steps are 
relatively wide. 

15 A tricritical point appears for an attractive interaction and bound states for a repulsive one. 
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facet width is proportional to A/vf- In Section [Vj it will be shown that the facet width is 
inverse proprotional to the effective correlation length on the surface: oc / S /A ~ A/vp- 
The interface profile near the edges of a facet will have a universal shape. Using Eqs. (^) 
and flUf ) one obtains: 

2 /int. 



/(o) 

» 

where the k° z = /(0)/A and = / S /A notations were used. That is 



K = ± k° x \l + ^h\k x )] , (17) 



k z — k° z — — ——== \k x kl\ 3 / 2 . (18) 

This equation establishes a universal (non-analytic) jump of the curvature, well-known in 
surface physics. The 3/2 exponent is the Pokrovsky-Talapov (1978) exponent while the tran- 
sition is called the Gruber-Mullins (1967) or the Pokrovsky-Talapov transition. The above 
description of the Pokrovsky-Talapov transition is mainly based on an effective nearest- 
neighbour model, however, as previously demonstrated, the Pokrovsky-Talapov exponent 
and correspondingly the transition is not destroyed by short-range interactions and not even 
by long-range interactions. The effects of the long-range interactions are analyzed in more 
detail the following Section, where it will be seen that the roughening transition remains 
unchanged, only a new phase may appear. 

On the rough side of the surface, however, a simple Taylor expansion cannot be per- 
formed. As presented in the Appendix; on the rough side the free energy will always have 
the form: 

f[h(k x )} - f(0) oc h 2 (k x ) + ... . (19) 

similar to Eq. ([T2]). 

Other characteristics of the roughening transition can be re-written for the Fermi surface, 
e.g., the universal jump in the curvature of the Fermi surface before and after the transition, 
or it can be shown that near the edge of the facets on the rough curved surfaces the steps 
are spontaneously generated with a square root density. 
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As presented in the Appendix, a Pokrovsky-Talapov transition appears by varying the 
a doping, or external magnetic or electric fields. As a function of temperature, however, 
such a transition is a Kosterlitz-Thouless transition if the width of the facet vanishes at the 
transition, otherwise is a standard second-order phase transition. 

VI. FERMI AND NON-FERMI LIQUIDS 

The conclusions from the previous Sections were that the rough phase of the Fermi surface 
corresponds to the metallic (Fermi liquid) state of a system, while a faceting (Pokrovsky- 
Talapov type) transition will lead to an instability of the Fermi surface corresponding to 
the appearance of some kind of condensate. In both cases, the surface may be characterized 
by the density (occupation probability) distribution p(k) oc 5n(k), where 5n(k) was defined 
in Eq. (0). As mentioned in Section ||, in the rough phase the density-density correlation 
function has an algebraic decay but an exponential decay for the smooth faceted phase. It is 
important to emphasize that these correlation functions appear in both translational and ori- 
entational order of the surface, which are decoupled (Haldane 1994), with the understanding 
that the orientational fluctuations cannot change the shape of the surface. 

The disordered (rough) phase is a non-interacting or Tomonaga-Luttinger liquid, that is 

. . . K , cos[7r k T h(k r )] , , 

WW = TT2 " K ,7 > (2°) 

7r k£ ki 

where K is the well-known roughness constant. Equation ( p0| ) can be viewed as the local 
slope operator so, the height-height correlation function 

([h{k x ) - h{0)) 2 } = 2 f X dk' x (k x - k y ) (p(k x )p(0)) , (21) 

J 

will have the well-known form: 

<[MM - MO)] 2 ) « *-£ in*. + |j , (22) 

with the universal constant K taking the K = 1 value, corresponding to the notations used 
in Eqs. (|20|) and (pip . It can be seen from Eq. ( p0|) that the non-interacting part adequately 
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describes the properties of surface. The second term of Eq. (|20|) which is the characteristic 
contribution of the Tomonaga-Luttinger models, is not affecting the height-height correlation 
function. Accordingly, the low-energy and long-wave length properties of the metallic phase 
are non-interacting in nature. By this we reconfirmed the Fermi liquid picture. 
The faceted phase on the other hand, maps into a Luther-Emery phase: 



where £ oc X/f s ~ up/A is the effective correlation length on the surface and A the quasi- 
particle gap. The / s is step free energy defined in Section [TV]. In Section [Vj, the above result 
was previously derived from a different point of view, namely using surface minimization 
arguments and determining the Lagrange multiplier of the surface. 

In the orientational order, i.e. in the k z direction, the situation is not so simple and 
detailed microscopic calculations are needed (Gulacsi 1996). However, we know that in the 
metallic phase (rough surface) there must be a step singularity^ in n(k) oc p(k) of the 
n(k) = 8(k — kp) form with 6(x) being the step function. While, in the smooth phase of the 
Fermi surface there are no singularities, e.g. in a mean-field description n(k) ~ l/[A 2 +v F k 2 ]. 

Concerning a non-Fermi liquid type behaviour, it seems (based on the above presented 
analysis) that there are no other phases left. The exact solutions of nearest-neighbour in- 
teracting lattice models, or surface models, exhibits three phases (see, Appendix). First, 
one which is ferroelectric, corresponding to a completely locked (or frozen) solid surface, 
corresponding to a ferromagnetic type transition, i.e., a Pomeranchuk (1958) type insta- 
bility of the Fermi surface^. The second, disordered phase and the third, ordered phase 

16 The presence of singularities as we cross a surface is not a surprise. In the case of a standard solid- 
liquid (or solid-gas) interface the mass density has a singularity at the surface bp = p &0 \\& ~ /^liquid; 
which determines the Lagrange multiplier defined in Section [n] (Nozieres 1991). 

17 The frozen ferroelectric phases of the vertex models are the only ones that exhibit true long-range 
order corresponding to a ferromagnetic state. 




(23) 
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corresponds to the rough and flat (faceted) surfaces, respectively. In this context, one of 
the most appealing possibilities to generate a liquid state, that is a metallic phase, but not 
a Fermi liquid, would be a phase intermediate between the rough and flat surfaces, which 
shares the properties of both. That is, a new phase between these two such that the new 
surface would be characterized by an exponential decay of the translational order but an 
algebraic decay of the correlation function of the orientational order. The existence of such 
a surface (sometimes called a hexatic fluid) is still controversial in the surface community 
(Bauer 1987), however, it would correspond to the disordered flat surface introduced by den 
Nijs (1989) which is in a pre-roughening state, corresponding to the presence of a "Haldane 
gap" type gap. Being so, it has what we could call hidden topological order. 

The appearance of such a phase is due to short-range interactions (Nijs, den 1989). 
In Section [V], I have shown that next-nearest-neighbour or short-range interactions cannot 
disturb or change the Pokrovsky-Talapov melting transition. However, what they can do, is 
to generate the above introduced disordered flat surface (Nijs, den 1989). 

An algebraic decay of the orientational order, means that n(k) will have some singularity 
as we cross kp, and the surface would correspond to a metallic phase. However, the metallic 
phase is not a Fermi liquid because the translational order is gapped. The step free energy is 
zero in such a non- Fermi liquid phase (Nijs, den 1989), so based on the definitions introduced 



in Section |TTT| the surface is "rough", but the height-height correlations are not divergent 
and thus the reason why such a phase is refered to as pre-roughening phase. 

Concerning the translational order, which is gapped, the surface exhibits flat surface 
facets similar to the previously analyzed faceted phase. The main difference is that, in the 



later phase the facets have a well-defined macroscopical width (see, Section |HI|) , while in a 
pre-rough phase only the average surface is flat. The surface structure is characterized by 
perfect antiferromagnetic order of the heights h(k x ) at random positions k x . Clearly, the 
surface structure is complex; that is why it is so difficult to describe analytically a non-Fermi 
liquid in higher dimensions. Because of the random positions of k x , no nesting can occur 
on the Fermi surface, accordingly no long-range order, or condensate either, and indeed, 
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the system will be metallic. The step free energy being zero, Eq. (|12|) is still valid, with a 
specific /(k,k'). Within the present phenomenological approach the exact form of /(k,k') 
cannot be determined. However, some insights can be gained following a simple analysis: for 
a Fermi liquid, as mentioned previously, both orientational and transversal correlations are 
algebraic. Being so, there is a singularity in n(k) and Z = {1 — <9XR e (k, u)[u = e(k)]/cta;} _1 
is finite. Let us consider the pre-rough, non-Fermi liquid, phase. As previously discussed, 
the orientational order is still algebraic; some sort of singularity will appear in n(k) as we 
cross hp. The translational order being gapped Z is zero. So, the emerging non-Fermi liquid 
will have a singularity similar to a Tomonaga-Luttinger liquid in one-dimension. However, 
I must emphasize that the above analysis is approximate within a semi-classical approach 
to the Fermi surface, and a microscopic treatment is needed (Gulacsi 1996) to prove exactly 
the above statement. 

VII. CONCLUSIONS 

Based mainly on a semi-classical approach, I have demonstrated that there is a one-to- 
one correspondence between the properties of a two-dimensional interacting electron system 
and its one-dimensional Fermi surface. Though, some of the evidence can be generalized to 
a three-dimensional interacting electron system. 

In a generic scenario, the Fermi surface can have three distinctive phases, namely the 
locked solid, which would correspond to a ferromagnetic type instability, a rough phase, 
which I have proven to be the metallic or Fermi liquid phase and a flat faceted phase, which 
is an instability of the Fermi surface against the appearance of a quasiparticle gap. By 
varying the interaction strength, external magnetic or electric, fields, or doping, we can 
cross these phases in the previously presented order, that is, frozen solid rough, melted 
surface •<=>- flat, faceted surface. It is interesting to note, that there is no transition between 
frozen solid and flat faceted surface. The transitions between these phases are: i) first or 
second-order from the frozen solid to rough phase, ii) Kosterlitz-Thouless or Pokrovsky- 
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Talapov melting between the flat faceted and rough surfaces. The properties of the rough 
and faceted surfaces and the transition between these two phase were described in great 
detail in Sections [Tl] - [TX|. In some cases, a third phase can appear between the rough and 
flat surfaces, namely a pre-rough surface with hidden topological order[^. So, the possible 
phases in the order of their appearance are: 

Locked •<=>- Rough •<=>- Pre — rough -<=>- FlatFaceted , (24) 

and these will correspond to 

Ferromagnetic •<=>- Fermi Liquid -<=>- Non — Fermi Liquid -<=>- Condensate . (25) 

The corresponding transitions are: i) first or second-order from the frozen solid to rough 
phase, ii) Kosterlitz-Thouless or Pokrovsky- Talapov melting between the rough and pre- 
rough surfaces. Hi) first or second-order between the flat faceted and pre-rough phases. 
Such a scenario predicts for two-dimensional models (e.g., the Hubbard model) a non- Fermi 
liquid behavior close to half filling, that is between the low-density Fermi liquid behaviour 
and the half-filled band case. 

In regards to the rough phase, a true analysis of its properties requires a discussion 
of its dynamics. In the Fermi liquid theory there are the sound modes and also solutions 
with imaginary frequency that represents the particle-hole continuum. Immediately, it can 
be seen that, the zero sound corresponds to the bosonic excitations of the one-dimensional 
system. In a one-dimensional system, we would not expect a continuum to appear because of 
reduced phase space. However, treating the Fermi surface as an effective equilibrium surface 
we can attach separate orientational and translational degrees of freedom, as demonstrated 
in the previous Section. Such that the orientational bosonic fluctuations are the sound 



18 In connection with the following discussion: the pre-rough phase will have sound modes (the 
orientational order remains unchanged compared to the rough phase), however, the particle-hole 
continuum will disappear because the translational order is different. 
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modes, while the translational fluctuations can generate a particle-hole continuum. Being 

a rough phase, these fluctuations are Gaussian, in a discrete systemp 9 ! the Lagrangian is: 

(g/4ir) J | V0| 2 , with (p = tt£/2 and i = ±1. We can identify this with the six-vertex 

model (see, Appendix) in the following way: in the eight-vertex model (Baxter 1990) the 

magnetic charge attached to the vertices seven and eight is m = ±1, due to which the scaling 

dimension of a vortex operator is gm 2 /2 = g/2 and, accordingly, the scaling dimension of 

its corresponding fugacity is 2 — gm 2 /2 = 2 — g/2. If the Boltzmann weights of vertices 

seven and eight are vanishingly small, the free energy is singular with an exponent inverse 

proportional to the fugacity which identified with the singular behaviour obtained by Baxter 

(Baxter 1990) gives : cos(g7r/4) = — cos /i, where cos/x is the anisotropy parameter of the 

model (see, Appendix). Thus, the translation fluctuations of the Fermi surface will have a 

ground-state corresponding to that of the planar six- vertex model (Baxter 1990): 

cos/i sinu f+°° , sinh(7r — a)x 

lu = " / dx — — ^ V — > ( 26 ) 

4 2 J-oo sinh tcx cosh fix 

with low lying states forming a particle-hole continuum bounded by 



7r sin jj, | , sin [i 



sin k I < lj < n 



k 

sin — 

2 



(27) 



2 fj, ~ ~ n 

Thus, not only bosons corresponding to zero sound will appear on the Fermi surface but 
also a particle-hole continuum. 
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19 As already presented in Section [T^, at large length scales the discreteness of the step heights 
is irrelevant, due to which the surface always can be described by a Gaussian model, having the 
properties of a transverse vibrating elastic network. 
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IX. APPENDIX: THE POKROVSKY-TALAPOV TRANSITION 



In the past twenty years, incommensurate periodicities have been observed in many phys- 
ical systems and have become the subject of considerable theoretical interest. Two incom- 
mensurate structures in contact can overcome their rigidity and form a common periodicity. 
This transition is known as the commensurate-incommensurate, or lock-in transition. 

Since its discovery in the context with dislocation theory (Frank and Merwe, van der 
1949), the commensurate- incommensurate transition had a major role (for a review, see Bak 
1982; or Nijs, den 1988) in understanding the physics of charge-density-waves, incommensu- 
rate lattice problems, the general problem of incommensurate periodicities, two-dimensional 
(XY) models, theory of phase transitions in magnetic systems and ID quantum systems, 
and the correlation driven metal insulator (Mott) transition (Gulacsi and Bedell 1994b). 

In most of the systems and models mentioned above, the problem of characterizing the 
commensurate-incommensurate transition has been reduced to the solution of the ID sine- 
Gordon model. The level of incommensurability is measured by a soliton density, n so i.(//), 
where /i is the chemical potential. The commensurate-incommensurate transition happens 
at /i — E so i, where E so \_ is the single soliton energy. For ji < E so i and n so i. = the 
system is locked in a commensurate phase, while for \x > E so i, i.e., n so i. ^ the structure is 
incommensurate. 

In the classical limit the commensurate-incommensurate transition occurs with the emer- 
gence of a soliton lattice in the soliton gap (Frank and Merwe, van der 1949; and Horowitz, 
1982). At low soliton density, the solitons repel each other, with an exponential type 
interaction (Horowitz 1982), and the transition takes place with a logarithmic behavior 
n so i. oc ln _1 (yU — -Esoi.), i-e., corresponding to a critical exponent (3 = 0. In the full quan- 
tum problem, it is known (Haldane 1982), that the soliton lattice will melt, giving rise to 
a Luttinger liquid which emerges at the top or bottom of the soliton gap. In this case the 
solitons will repel each other much stronger than in the classical limit, with an interaction 
proportional to n^ ol , in the low soliton density limit. The incommensurability close to the 
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transition behaves as n so i, oc (// — A) 1 / 2 , where A is the soliton gap, corresponding to a 
critical exponent (3 = 1/2. This is the Pokrovsky-Talapov exponent (Pokrovsky and Ta- 
lapov 1978). In some applications this transition is characterized with the variation of the 
thermodynamic potential, f2 = E(n so i) — Lin so i, for which we obtain cx rig 0l = (/x — A) 3//2 . 
From the experimental point of view, the relevant quantity to observe this transition is the 
charge susceptibility, Xsoi. = (9///9n so i.) _1 , which is divergent at the transition Xsoi. oc l/"n. so i.- 
This divergent behavior implies a divergence of the effective mass of the charge carriers, and 
accordingly the divergence of the specific heat coefficient. Thus, the basic characteristics of 
the Pokrovsky-Talapov transition are: 

n sol .(fx) cx (fx - A) 1 / 2 , Q(fj,) cx (/i - A) 3 / 2 , Xsol .(/i) cx (li - A)' 1 ' 2 . (28) 

Or, as a function of the soliton density (incommensurability) are: 

//(n so i.) cx A + n 2 ol , E(n soh ) cx An so i. + "L./ 2 > x(n so i.) oc l/n so i. . (29) 

The Pokrovsky-Talapov transition, however, is typical for all integrable models. Its im- 
portance in connection to the 2D XY, Gaussian, Coulomb-Gas models but, most important 
to the equivalent (Chui and Weeks 1976; Beijeren, van, 1977; and Jayaprakash and Saam 
1984) solid-on-solid surface models led a to thorough analysis of its properties. Recognizing 
(Nijs, den 1981; and B. Horowitz, Bohr, Kosterlitz and Schulz 1983) the equivalence be- 
tween the Gaussian models and the direct field six-vertex model (defined below) led to the 
remarkable result that the six-vertex model in a direct field is an ideal model to describe the 
commensurate-incommensurate transition. The direct field in the six-vertex model corre- 
spond to the chemical potential in Eqs. ( p8[) and (p9|), while the polarization of the six- vertex 
model is equivalent to the soliton density, n so i.. 

The five- and six- vertex models are the simplest, exactly solvable model describing surface 
transitions. In general, the vertex models define a statistical mechanics on a (for simplicity) 
square lattice as follows. On each link there is a degree of freedom taking two values which 
is represented (traditionally) by an arrow, i.e. — > and <— . To each vector therefore there 
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corresponds 16 possible configurations. To each of these we associate a Boltzmann weight 
Ui = exp(— Ei/ksT). This general model is referred to as the sixteen-vertex model. The 
exact solution of this general ferroelectric model is not known. If we choose the Boltzmann 
weights, so that six vertices with two entering and two exiting arrows (or eight vertices 
including also the cases where four or zero arrows exit) have finite weight, then we define 
the six- or eight-vertex models. The six-vertex model is exactly solvable via the algebraic 
Bethe Ansatz, for details see Lieb and Wu (Lieb and Wu 1972) and Baxter (Baxter 1990). 
The exact solution has at most four independent variables (this is because the zero energy 
can be chosen arbitrary and the ice rule implies uj^ = uj^)f^ In this situation the model 
is called the asymmetric six-vertex model, in contrast to the symmetric six-vertex model, 
where the number of independent variables are two. The symmetric model in the presence of 
h and v external fields is equivalent to the asymmetric one. The energies of the six possible 
arrow configurations of the symmetric model are E\ = £2, £3 = £4 and £5 = e§. Similarly, 
the field-free five- vertex model is defined by three independent Boltzmann weights, e.g., e^, 
62 = 00, £3 = £4 and £5 = Eq and its exact solution can be determined via the algebraic 
Bethe Ansatz (Gulacsi, Beijeren, van and Levi 1993). 

1.) For the six- vertex model, the exact free energy obtained with the algebraic Bethe 
Ansatz (Lieb and Wu 1972; and Baxter 1990) exhibits three different analytic forms for 
the three different phases, i.e., two frozen-in ferroelectric phases (with the polarization, 
p = ±1), an antiferroelectric phase (p = 0) and a disordered phase (paraelectric, p = 0). 
The transitions to the ferroelectric phases is of first-order, while the transition from the 
paraelectric phase to the antiferroelectric one is of Kosterlitz-Thouless type. 

Introducing a vertical field, v in the symmetric model is still easily solvable with the 
algebraic Bethe Ansatz. All the above mentioned transitions are second-order now, as a 
function of temperature and Pokrovsky-Talapov as function of the external field. The largest 



20 We use the traditional ordering of the vertices (Baxter 1990). 
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eigenvalue of the transfer matrix is determined from the Bethe Ansatz, z(p) in notation of 
(Lieb and Wu). The exact free energy is obtained as —F/k B T = max z{p) -i< p <i- All the 
three transitions are similar in nature, accordingly we describe only the antiferroelectric one, 
p —>■ 0. The free energy in the ordered phase (smooth faceted surface) is determined by (Eq. 
(304) from Lieb and Wu 1972): 

z{p) oc z'(0)p + lz"'(0)p 3 + ... , (30) 
o 



with z" = 0, i.e., there is no second-order contribution. Eq. ( pOf) is identical to Eq. (p9|) . 
Close to the transition p oc (v - A; B T|^(0)|) 1/2 (Eq. (340) from Lieb and Wu 1972), tt oc 
(v — kBT\z'(0)\) 3 / 2 (Eq. (342) from Lieb and Wu 1972) and the magnetic susceptibility 
X oc (v — ksT\z' '(O)!) -1 / 2 which is equivalent to Eq. (p8|). The same result is obtained for 
the full asymmetric six-vertex model, i.e., with both h and v being present, however, due to 
the complexity of the resulting equations numerical analysis is required. 

In the disordered phase (rough surface), however, the free energy has the form (Eq. (316) 
from Lieb and Wu 1972): 

z\P) « — ^— cos (^)P + ••• ' ( 31 ) 
where cos/i is the anisotropy parameter of the model and = (1 + u>ie lfi ) / (e ltl + uji). 



In connection with Section |V| it is important to emphasize the fact that no linear term 
appears in Eq. (^). In this part of the phase diagram (i.e., on the rough surface) Taylor 
expansions cannot be performed, as compared to Eq. (|30|), result of the Eq. (|31"D or (|33D 
type are obtained directly from an integral equation. 

2.) In the five-vertex model, in contrast to the symmetric six-vertex model, the arrow- 
reversal symmetry is broken. Again there are three different analytic forms of the exact free 
energy (Gulacsi, Beijeren, van and Levi 1993) corresponding to two frozen-in ferroelectric 
phases (p = ±1), one antiferroelectric phase (p = 0) (similar to the symmetric six- vertex 
model) and a ferrielectric (—1 < p < 1) phase. All the transitions are of second-order. From 
all these phases, the ferrielectric and ferroelectric phases are connected to surfaces (see later) 
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due to which I briefly present their critical behaviour. In this case, e.g., p — > 1 and the free 
energy in the ordered phase (smooth facets) is (Eq. (A7) and (A8) from Gulacsi, Levi and 
Tosatti 1994): 



/(1-p) oc f(0) + k B T 



-^ln(^U21n^ 3 



P 



2 



_ £ g4j + ^) h ,^ + _ _ (32) 



6 (u>i — LU3) 3 \ 2 

As in the six-vertex model, see Eq. (|30| ) there is no second-order term. Close to the transition 
1 — p, Q and the susceptibility behaves similar to the previously presented six- vertex case 
or Eq. (PSD- The free energy on the ferrielectric side (rough surface) has a particular simple 
form in the non- interacting, zero external field limit: 

/(1-p) oc /(0) + ^k B T(^) 2 + ... , (33) 

while for the full formula see Eq. (21) from Gulacsi et al. (Gulacsi, Levi and Tosatti 1994). 
Because of the simplicity of Eq. ( |33"D the five-vertex model is most adequate to model the 
Pokrovsky-Talapov transition. 

As a third example, I present briefly a variant of the Pokrovsky-Talapov transition ap- 
pearing on surfaces, similar to Section [V], known also as Gruber-Mullins (1967) transition. 
Surfaces of simple cubic crystals are described by the five-vertex model (Gulacsi, Levi and 
Tosatti 1994), while surfaces of body centered crystals with the six- vertex model (Beijeren, 
van, 1977; and Jayaprakash and Saam 1984) As presented in Section [T|, knowing the free 
energy of the surface we can determine the Cartesian coordinates, z(x, y), of the correspond- 
ing surface from the Wulff construction. The free energy of the surface we obtain directly 
from the exact solution of the five- or six-vertex model, depending on the crystal structure. 
For the simple cubic surfaces, p — > 1 and 1 — p oc (x — xq) 1 ^ 2 , while for surface of body 
centered crystals p — > and p oc (x — Xq) 1 ^ 2 . In both cases, z — z oc (x — Xo) 3 ^ 2 and the 
average spacing between steps oc (x — xo) ~ 1 ^ 2 , equations identical to Eq. (|28|). 

As mentioned above, a general characteristic of the (Luttinger) liquid phase preceding 
the Pokrovsky-Talapov transition is that its quasiparticles strongly repel each other. The 
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basic equation obtained for the ground-state energy in Eq. (p9| ) can be interpreted as follows: 
the linear term in n so \ sums up the energy contributions of the non-interacting quasipar- 
ticles (solitons, holons, surface steps, etc., depending on the model under consideration). 
The missing quadratic term in n so i. shows that there is no effective attraction between the 
quasiparticles (see below) and the n^ ol term suggests of the presence of strong repulsion. 
This repulsion is discussed in detail in Section [V|. 

As all the previous proof were given for an effective nearest-neighbour interacting sys- 
tems, it is very important to realize that the Pokrovsky-Talapov transition is very stable 
with respect to different types of interactions. Short range interactions, with range much 
less than the average distance between the quasiparticles contribute to the ground-state en- 
ergy by an order of n^ ol thus, leaving the Pokrovsky-Talapov transition unaltered as it was 
proven in Section |V|, where also the cases of long-range interactions was analyzed. 

By the above examples we want to emphasize that the integrable models have common 
ground-state properties^] thus, a Pokrovsky-Talapov type transition will exist in all related 
models. 



It should be emphasized that the above mentioned mappings between different models estab- 
lishes an equivalence between the eigenvalues of the different models. That is, the ground-state 
energies and thus, the thermodynamic properties of the different models can be made equivalent. 
However, non of these mapping schemes preserves the correlation functions. 
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XII. FIGURE CAPTION 

Fig. I. The surface of a square (solid line) and its corresponding Wulff plot (dashes 
line). 
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Fig. 1 



